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Abstract. It has been recently discovered that exceptional families of Sturm- 
Liouville orthogonal polynomials exist, that generalize in some sense the clas- 
sical polynomials of Hermite, Laguerre and Jacobi. In this paper we show 
how new families of exceptional orthogonal polynomials can be constructed by 
means of multiple-step algebraic Darboux transformations. The construction 
is illustrated with an example of a 2-step Darboux transformation of the clas- 
sical Laguerre polynomials, which gives rise to a new orthogonal polynomial 
system indexed by two integer parameters. For particular values of these pa- 
rameters, the classical Laguerre and the type II X^-Laguerre polynomials are 
recovered. 



1. Introduction 

It has recently been shown that the hypotheses of Bochner's Theorem [2] on 
the characterization of orthogonal polynomial systems defined by Sturm-Liouville 
problems [17] can be relaxed to give rise to new complete orthogonal polynomial 
systems. In some recent papers, the concept of an exceptional polynomial subspace 
[71 [S] and the closely related notion of exceptional orthogonal polynomials 112] 
have been introduced. Like the classical examples, exceptional orthogonal polyno- 
mials are eigenfunctions of a second-order differential operator, but the sequence 
of eigenfunctions need not contain polynomials of all degrees, even though the full 
set of eigenfunctions still forms a basis of the weighted L^ space. 

The first examples of such polynomials were the Xi-Jacobi and Xi-Laguerre 
polynomials introduced in [TTJ [T^] . In this co-dimension 1 full character- 

ization of all Sturm-Liouville polynomial systems is available thanks to the clas- 
sification of codimension one exceptional polynomial subspaces performed in [11) . 
The first relation between exceptional polynomials and the Darboux transformation 
was made by Quesne in the context of shape invariant potentials in supersymmet- 
ric quantum mechanics [191 120] . where examples of X2 polynomials are also found. 
Higher codimension families were described for the first time by Sasaki and Odake 
[2T] for arbitrary codimension £. These authors also realized [22] that two distinct 
families of X^-Laguerre and X^- Jacobi polynomials exist, now labelled as type I and 
type II. The existence of precisely two families for each class was explained in [T3] 
for the Laguerre case, by showing that these families are obtained from the classical 
ones by means of an isospectral algebraic Darboux transformation |9J . Indeed, there 
are four families of algebraic factorizations, only two of which give rise to isospec- 
tral transformations [10]. The equivalent result for the Jacobi case was shown in 
[14] . and similar results on the relation between the Darboux transformation and 
exceptional polynomials have been published in [24] . 
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Exceptional orthogonal polynomials and their associated exactly solvable poten- 
tials have found numerous recent applications in a range of problems in mathe- 
matical physics, such as mass-dependent potentials [TH], Fokker-Planck and Dirac 
equations |l6j, supersymmetric quantum mechanics (5] [15] or quasi-exact solvability 

m- 

A natural question that arises from the previous works is whether the isospectral 
Darboux transformation can be iterated in a Darboux-Crum form 4, 3 in order 
to obtain families of exceptional orthogonal polynomials that enjoy spectral char- 
acterizations and completeness properties similar to the ones described above. Our 
purpose in this paper is to give a positive answer to this question by showing that 
two-step isospectral Darboux-Crum transformations can be used to construct a new 
set of exceptional orthogonal polynomials of Laguerre type. This family is indexed 
by two integer parameters (mi,m2) and a real parameter k. The sequence starts 
with degree £ = mi + m2 — 1 and it is complete in a weighted space endowed 
with a non-singular weight. Particular choices of mi and m2 correspond to the 
type II Xf-LsLgiierie and to the classical Laguerre polynomials. 

Our paper is organized as follows. In Section[2l we define the two-step exceptional 
Laguerre polynomials, present their characterization as the unique polynomial so- 
lutions of a second-order linear differential equation (Theorem [T]) and state their 
completeness properties (Theorem[2|). We devote Section[3]to the proof of Theorem 
[1] using rational factorizations, an isospectral Darboux-Crum transformation and 
the higher-order intertwining relation (|37p . Section |4] addresses then the complete- 
ness property stated in Theorem [2] through a polynomial approximation argument 
for weighted L^ spaces. 

2. TwO-STEP EXCEPTIONAL LAGUERRE POLYNOMIALS 

Our goal in this section is to introduce the two-step exceptional Laguerre polyno- 
mials, present their relation to second-order linear differential operators and state 
their completeness and orthogonality properties. The proofs of these results will be 
given in Sections [3] and H] 

We begin by recalling that for any given real number fc, the classical nth degree 
Laguerre polynomial Ln{z) = l\i\z) is the unique polynomial solution of the 
second-order differential equation 

(1) Ck[Ln] = -nLn, 
where 

(2) Ck[y]=zy" + {k + l-z)y', 
normalized by the condition 

(-z)" 

Ln(z) = ; h lower degree terms. 

nl 

For fc > — 1, the sequence Q spans a dense subspace of the Hilbert space 

(3) L2([0, oo),Wkdz), Wkiz) := z'=e-^ 

and forms an orthogonal polynomial family satisfying 

roo ( r(fc+«+i) . 

(4) / L^)Lfw,{z)dz=\ 
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The two-step exceptional Laguerre polynomials are now defined in the following 
way. Let k be real and < mi < m2 integers. Set 

(5) £ — mi + m2 — 1 
and define 

(6) r?i2 = W[r?i, 772], where = ii7„''~^\ a =1,2, 

and where W denotes the Wronskian operator. Consider now the second-order 
linear differential operator C defined by 

(7) Ck[y] = zy" + (fc + 1 - z)y' + —{r^'i^iy ~ y') + WK, 772]^)- 
Let n > ^ be an integer. Set 

(8) C := ((toi - m2){k + 2 - mi + n)(fc + 2 - m,2 + n)y^ 
and let L„ be the polynomial defined by 

(9) L„ := ^'=.'"1^'"=) = Cz-''W[m,m,^^^^L^n-I\ « > ^ 

(-z)" 

(10) = 7 -rr-. ; 7 + lower degree terms. 

(n — iy. mil 7712' 

For a lighter notation, we shall use L„ instead of ^^'^^'"i'™^) throughout the paper. 
Whenever the symbol I/„ is used on its own, the rest of the parameters G R, 
TOi, ?Ti2 G N are understood. 

The main results of our paper are the following. 

Theorem 1. The above-defined i„ is the unique polynomial solution of the differ- 
ential equation 

(11) Ck[Ln] ^ {2i - n)Ln, n>i 
subject to the normalization (jlOp . 

Theorem 2. If k > m2 — 2, then 7712(2) > for z > 0. Furthermore, the above 
defined polynomials {^nlj^^ span a dense subspace of the Hilbert space 

(12) nk,Tni.Tn2 ■^'i^'^{[0,Co),Wk,mi,m2dz), #fc,„ii (z) := Z*" 77^2^6"^ 

and satisfy the following orthogonality relations: 

CT{k + n + 3~i) 



(13) / LnLjWk,mi,m2dz = <^ {mi - m2) {n - ey. 

If ?7ii = 0,777,2 = 1, then by well-known identities for the Laguerre polynomials 

f (.,0,1), , _ z^Lj^^'^'iz) + (fc + 2)(2zLi'-+^)-(z) + (fc + l)Lit+'\z)) 
" ^ ^" (fc + 2-f 77)(fc-M + n) ~ " 

Therefore, the above-defined polynomials generalize the classical ones. If 77ii = 
0, 7712 = 777-1-1, where 777 > 1, then 

-—^ L('=."-+i) = z-^^W[47^-i), z'^+14'^+i)] = Lt,.„, n > m, 

k -\- Z -\- n — m 

where the latter are the recently introduced exceptional Laguerre polynomials of 
type II mi [13] . Therefore, the newly defined polynomials L„ also generalize and 
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extend the previously described class of exceptional Laguerre polynomials. This is 
natural since for mi — one of the factorizing functions becomes a constant and 
therefore we are dealing with a 1-step Darboux transformation. The two families 
of Xf-Laguerre exceptional polynomials of arbitrary codimension were shown to be 
obtainable by means of a 1-step Darboux transformation from the classical ones 

3. Rational factorizations 

This Section is devoted to the proof of Theorem [I] characterizing the two-step 
exceptional Laguerre polynomials L„ as the unique normalized polynomial solutions 
of the differential equation pD) . 

Let 

T\y\ = v{z)y" + (l{z)v' + r{z)v 
be a second-order differential operator with rational coefficients qiz) and r(^). 

Definition 1. A rational factorization T\y\ is an operator identity of the form 

T ^ BA+~X 

where A[y],B[y] are first-order differential operators with rational coefficients. 
We call 

f = AB+~X 

the partner operator corresponding to the above rational factorization. Note that 
the following formal intertwining relations therefore hold by construction: 

(14) AT = TA, TB = BT. 

Given a rational factorization, let (j){z) be the unique (up to constant multiple) 
solution of A[(f)] — 0. Observe that w{z) := (p' (z) / (j){z) is rational and that T[(j)] = 
\<j). For this reason we will call (^(2) a quasi-rational factorization function. Also 
note that the following identities hold for A and B: 

(15) A[y]^^-^W[4>M^b{z){y' -w{z)y) 

(16) B[y] ^^S^W[ly]=h{z){y' ^w{z)y), 

where h{z) is a rational function called the factorization gauge [T3|, and 

I P - "1 Z ( f 

b——, w ~ —w 1 , = exp w dz 

b p b \J 

The partner operator has the form 

(17) f[y]^piz)y" + qiz)y' + f{z)y 
where 

(18) q Q + P f ^ -p{w + li'^) — qw + X. 

Let us also define measures Wdz and Wdz by setting 

W =p-'^exp(^J ^dz^ , W = ^W. 
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The above defined operators A and —B are formally adjoint of each other with 
respect to these measures, in the sense that: 

(19) J A[f]g Wdz + J B[g]f Wdz = bWfg. 

It is important to note that a rational factorization is just special case of the well- 
known Darboux transformation for Schrodinger operators, albeit expressed relative 
to a general coordinate and gauge. Indeed, the operator T[(f>\ can be related to a 
Schrodinger operator 

(20) H[^] = + U{x)^P 
by means of a change of variable z = C{x) satisfying 

(21) Pia^)) = ci^f 

and a gauge transformation 

(22) -(/.T[0])oC = iTM, 

(23) = ificj)) o C, 

(24) fi = exp / — — —dz, 

J 4p 

(25) C/--(-C) + ^ 

In the physical variable and gauge, the above factorization corresponds to 

H={d, + u)i-d, + u)-~X, uix)^^, 

and is related to the factorization function (l){z) by (|22p . The partner operator in 
the physical gauge is given by 

(26) H = {-d, + u){d^ + u) - ~X = -d^, + U 

where the partner potential is given by the well-known Crum potential formula 

(27) U = U -2d,,\og^. 

Just as the Darboux transformation can be iterated to obtain the more general 
Darboux-Crum transformation, rational factorizations can be iterated and give rise 
to higher order intertwining relations. Indeed, we shall see that the operator L[y] 
defined in Q is obtained from the classical Laguerre operator L[y] by means of 
a 2-step factorization. We describe this construction in more detail, and thereby 
furnish a proof of Proposition [T] 

Fix k, nil, 7712 as above and let 

(28) To = Ck+2 

be the ordinary Laguerre operator. The quasi-rational functions 

(1), -.^ z-''-^ri„ i^l,2 

are factorization functions of To with eigenvalues Xi = k + 2 — rm O Section 6.1]. 
Consider the rational factorization of To corresponding to the factorizing function 
4>i and factorizing gauge 

(29) bi = zr]i 
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which gives 

(30) To = BiAi + ~Xi, 

where by equations (fT5|) and (|16|) it is clear that 



(31) ^iM = ^>V[0i,y], BM^y^. 
Let 

(32) Ti := To = A^B^ + Ai 
be the corresponding partner operator. Note now that 

<^12 Ai[02] = 02] = z-^-'rii2 

is a quasi-rational factorization function of Ti with eigenvalue A2 = k + 2 — m2. This 
means that 0i2 can be used as factorization function for a rational factorization of 
Ti, together with the factorization gauge 

(33) bi2 = . 

We thus obtain the factorization 

(34) Ti =^2^2 + A2 
where 

A2[y]-^W[0i2,y]. 

012 

We call T2 the partner operator of this last factorization 

(35) T2 fi = A2S2 + A2 
and we observe that 

(36) iA2A,M = h2h ^J^f'/''f = z-'^w[m,m,z'^'y]. 

>V[01,02j 

By construction, the following higher-order intertwining relation holds 

(37) T2A2Ai = A2A1T0. 

From the definition of the L„ polynomials in ^ and the above expression l\'6G\i we 
see that L„ are eigenpolynomials of the differential operator T2: 

(38) T2[Ln]^{i~n)Ln. 

To conclude the proof of Theorem 1 , we must establish that T2 and Ck are related 

by 

(39) A - T2 + e. 

Let Qi and r^, i = 0, 1, 2 be the coefficients of operators Ti above: 

(40) T,[y] = zy" + q,y' + ny, ^ = 0,l,2. 
From ([2]) and ((28|) . we see for instance that 

(41) qo = k + 3-z, 

(42) ro = 
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and hence bv (fTSl) ([29l) (|33l) . 

(43) q2=qi+p' -2p^ 

= go + 2p'-2plog(6i6i2)' 
= l + fc-z-2z^. 
which matches the first order part of ([7|). It remains to demonstrate that 

(44) r2 = —{ri[2 + mv'i,V2])-e, 

and for this purpose it will be useful to compare these expressions in the Schrodinger 
gauge by making use of the Crum potential formula. 

Let Ui{x), iii(z), i = 0, 1, 2 be the potentials and gauge factors corresponding to 
the operators Ti, as given by (|24)) - (|25|) . In particular, by p3t . 

(45) fio^e-^z^+^, 

(46) fl2 = 

Z1112 Vl2 

where the latter follows by (|43p. Since p{z) — z, the required change of variables is 

(47) C(^) = ^, 
and by the Crum potential formula (|27p. we have 

(48) U2 = Uo- 29,,logW[V'i,V'2], 
where ipi = (mo^i) ° C- It follows that 

(49) logW[V'i>2] - {logW[(j)i,(j)2] + ilogp + 21og^o) oC 

= log(r/i2oC)-(l + fc)logC-C 
By (HSl), (gS]) and (1121) we have 

.2 o C = 2(log W[^„^2])" + ^P^(-) - ^f^i^)- 

[fi2 o Q (a<o o C) 

Applying the identity 

/"// - 97.9 - (log(//.g))" + (log(//5))'(log(/g))' 
with / = /i2 o C g = fiQ o and relation (j46p gives 
r2 o C = (log(77i2 o C) - 2C - (2fc + 3) log C)"+ 

+ (log C + log(r/i2 o C))'(log C + log(?7i2 o C) - 2 log(Mo o C))' 

(^?12 ° 0" 



- (2C + (2A; + 3) logC)" + (log C)'(log C - 2 log(/io ° C))' 
'7i2°C))'(logC-(A*o°C))' 
+ 2(log(77i2oC))'(logC-(Mo°C))' 



{Vl2 ° C) 

+ 2(log(77i2oC))'(logC-(A*o°C))' 
{Vl2 ° 0" 
{Vl2 O 
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Now switching back to the algebraic variable gives 

^ V12 \z Mo/ V12 

(50) =£!?31 + (^_fc)^. 

f?12 ??12 

Since y = rji, i — 1,2 are solutions of the Laguerre differential equation, 

(51) zy" — {1 + k + z)y' + rriiy ~ 0, 
we have 

zri'l^ - (/c + z)7y;2 - 2zW[r;;, 77^] + £7^12 = 0, 
which together with ((50)) implies p4)) . the required conclusion. We have thus es- 
tablished that the polynomials L„ defined by ^ are obtained from the associated 

ik') 

Laguerre polynomials L„ by use of a two-step rational factorization. The proof of 
Theorem 1 is achieved by showing that the differential equation Ck in ([7]) satisfied 
by the polynomials L„ is essentially the partner of Laguerre's operator, correspond- 
ing to the iterated Darboux transformation. It is important to stress that the two 
factorizing functions of the iterated rational factorization: 



H=z 



are indexed by two integers mi and m2 that can be freely chosen. This fact should 
be put in comparison with the restricted choice of factorizing functions in the 
original Crum method [3] or its modification by Adler [1]. The main difference lies 
also in the fact that in the former constructions the Darboux transformations are 
state-deleting, while in this paper the factorizing functions correspond to isospectral 
transformations. 



4. Proofs of completeness and orthogonality properties 

Lemma 1. Let m2 > mi > be integers and k > m2 — 2a real number. Then, 
^I2{z) > for z > 0. 

Proof. Since y — r]i, i — 1,2 satisfies (|5ip . 7712 satisfies the inhomogeneous first- 
order differential equation 

W12 - (fc + 1 + z)rii2 + (m2 - mi)r]i7]2 = 0. 

This gives 

/•oo 

?7i2(z) = (m2 - mi)e''z^+'' / e-H-^-''f]i{t)f]2{t)dt. 

J Z 

If fc > m - 2, then by Theorem 6.73 of [25], L^''"^^(z) > for z > 0. The desired 
conclusion follows immediately. 

Let V denote the vector space of univariate polynomials, and Vn the subspace 
of polynomials of degree < n. Let £ = mi + 777,2 — 1 and let f„ C Vn+e denote the 
span of Li, . . . , Ln+g. Let Zi, i = 1, . . . , £ denote the zeroes of 7712(2). 

Lemma 2. A polynomial y e Vn+g belongs to if and only if 

(52) - 2z,y'{z,) +(z,-k+ ''f^^'^'^A y{z,) ^ 0, i = !,...,£ 

\ Vl2(^i) / 
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Proof. Suppose that y e £n- By equations (|43|) ([50l) of the previous section, 



Res 



— yy 

??12 f?12 



0, 



These are precisely the relations (|52p . Since the codimension of f„ in P„+f is 
exactly the constraints ([52]) define the subspace f„, as claimed. 



Lemma 3. Let r]{z) be a polynomial such that r]{z) ^ for z > 0. If k > —1, then 
the subspace rjV := {r](z)p{z) \ p G V} is dense in Hk- 

Proof. Since V is dense in Hk [IS], it sufhces to prove that V is contained in the 
closure of ijV. Let p and e > be given. We seek a polynomial q £V such 
that \\p — qriWk < e. To this end let zo = 1 + max{|z| : 77(2) = 0} and set 



r(z) 



' p{z-zo) 





Z > Zq, 

otherwise 



Set I ^ k + 2 deg rj and observe that 

z'^\rj{z)\^ < {z + zoY, 
Since l/|7;(z)| is bounded for 2; > 0, we have 



p{z - zo) 



ri{z - Zq) 



z > 0. 



p(z) 



ri(z) 



[z + zo)'e ""dz < 00 



Since V is dense in He, there exists a. q G V such that \\r — q\\i < e. Hence, 
q{z) ~ q{z + Zq) is the desired polynomial because 



\\p-qm\k 



\r{z + zo) - q{z + zq)Y \tj{z)\' z^e'' dz 



< I |r(z + zo)-g(z + zo)r|z + zore ""dz 
10 



< e 



'° I \riz)~q{z)\^z^e-'dz 




< e 

This concludes the proof of Lemma [3l 

We now use Lemma[2]and Lemma[3]to prove Theorem[21 The positivity of ?7i2(^) 
was established in Lemma [T] It therefore follows that the measure Wk.m-i,m2{z)dz 
has finite moments. Let us now prove that the set of all 2-step Darboux exceptional 
Laguerre polynomials 

£ = y = span{L„}J^^ 



is dense in 'Hk,mi,m2- Let / e 'Hk,mi,m2 and e > be given. Set / = //7?i2 and 
note that / G Hk because ||/||fe,mi,m2 — ll/llfc- Hence, by Lemma|3l there exists a 
polynomial p G V such that 



\\.f - PVuWk.rm, 



\f{z)-7^,2{z)p{z)\'z'e-'dz<e 



By Lemma m ^77^2 G which establishes the claim. 
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To conclude the proof of Theorem [2] let us prove the orthogonality (fT3|) of the 
exceptional polynomials. Let ^42/] 7^42/]; * — 1,2 be the first-order operators 
defined in the preceding section. For f,g eV, by (OITID ^ ^ 

{A2A,f,A2A,g), „^^ ,^^ - {B,B2AM,g),^, 
= (Si(Ti+A2)Ai/,g\ 

\ / k+2 

= ((To + A2)(ro + Ai)/,g) 

\ / fc+2 

The desired conclusions follow by ^ and ([4]). 

5. Summary and conclusions 

We have shown how the iterated Darboux or Darboux-Crum transformation can 
be used to construct new families of exceptional orthogonal polynomials, thus en- 
larging the class previously described in the literature. Due to the many properties 
in common with classical orthogonal polynomials, a systematic classification of all 
the exceptional families is an important goal, and the Darboux-Crum construction 
described in this paper seems an appropriate tool to achieve it. 

In the Schrodinger gauge and physical variable, this construction gives rise to 
new exactly solvable potentials in quantum mechanics. 
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